Nodes vs. minima in the energy gap of iron-pnictides from field-induced anisotropy 
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We develop the formalism for computing the oscillations of the specific heat and thermal transport 
under rotated magnetic field in multiband superconductors with anisotropic gap and apply it to 
iron-based materials. We show that these oscillations change sign at low temperatures and fields, 
which strongly influences the experimental conclusions about the gap structure. We find that recent 
measurements of the specific heat oscillations indicate that the iron-based superconductors possess 
an anisotropic gap with deep minima or nodes close to the line connecting electron and hole pockets. 
We predict the behavior of the thermal conductivity that will help distinguish between these cases. 
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Discovery of superconductivity in iron-pnictides re- 
energized the effort to understand the properties of the 
paired state in correlated electron systems. Materials 
containing FeAs [l|, Q or Fe(Se,S,Te) @, 3 layers are 
quasi-two-dimensional, and most are antiferromagnetic 
(AFM) at stoichiometry, which led to early comparisons 
to the cuprate superconductors (SCs) . Unlike cuprates, 
the pnictides remain metallic \MR in the AFM state, 
suggesting that the itinerant correlated picture is more 
appropriate. This view is also supported by the agree- 
ment between the band structure calculations and the 
measured Fermi surface [§-[l2j]. 

The complication is that all five d-orbitals of Fe ions 
contribute to the density of states (DOS) close to the 
Fermi level resulting in multiple Fermi surface (FS) 
sheets. The complete description of these materials in- 
cludes two to three hole bands at the center (r-point) 
of the Brillouin Zone (BZ) and two electron bands in 
the corner of the BZ, hereafter referred to as M-point. 
This multiband nature is essential for the ongoing de- 
bate about the structure of the superconducting order 
parameter in iron-based SCs. 

Superconductivity in pnictides is likely due to the mag- 
netically assisted electron scattering between the nearly- 
nested hole (h) and electron (e) FS sheets @, leading 
to a so-called state, with both pockets fully gapped, 
and A e = — A^. Detailed description of spin 

fluctuations and intra-band Coulomb scattering favors 
anisotropic gaps: only on the electron sheets for the A\ g 
(extended s-wave) representation, flfjl - 1201 ] , and on all FS 
sheets for the gap of B\ g (d-wave) character 2jJ, |22J . The 



latter gap shape is unlikely since ARPES measurements 
see nearly uniform gaps on hole FS. 1(1 23 - 28| 

The magnitude of the anisotropic component in the 
A\ g state depends on the values of the interaction pa- 
rameters and hence is material-dependent. Thus, possi- 
bilities range from isotropic A e , to a gap with deep min- 
ima on the electronic FS along T-M line, to a state with 



a pair of "accidental" nodes near this line [181 Il9l |29| . 
In the unlikely case of the dominant anisotropic compo- 
nent, the nodes move to positions along the sides of the 
crystallographic BZ. Signatures of low-energy excitations 



were found in Co-doped Ba(FeAs) 2 [3fJ, LaFePO[31 
BaFe 2 (Asi_ x P K ) 2 jHHsl and Fe(Se,Te) [I 



materials. 

However, the detailed gap structure of the pnictides, in- 
cluding the location of the possible nodes on the electron 
sheet, still needs to be unambiguously determined. 

Oscillations of the thermodynamic and transport coef- 
ficients in SCs with anisotropic gap as a function of the 
relative orientation of the magnetic field and the nodal 



(or quasi-nodal) directions [37H41I ] are extensively used 



to determine the position of the symmetry- enforced gap 
nodes [42|, |43| . The key prediction of the inversion of the 
anisotropy |39| . the switch from the minima to maxima 
for the field along the direction of the smallest gap, was 
recently confirmed |44l . A similar test was suggested for 
pnictides in Ref. [451 ]. and very recent measurements of 
the specific heat in the vortex state of Fe(Se,Te) j4(| were 
interpreted as leading to a surprising conclusion that the 
nodes of the gap are along the principal directions in the 
Brillouin Zone. This stimulated our study. 

In this Letter we develop the formalism for com- 
puting the properties of the vortex state of multiband 
two-dimensional (2D) superconductors under an in-plane 
magnetic field. We specifically address the states without 
symmetry-enforced nodes, such as pnictides. (4?J We an- 
alyze the behavior of the specific heat, (C), and the elec- 
tronic thermal conductivity, (k), focusing on the regime 
where the inversion of the oscillations occurs, and com- 
pare it with data in Ref. (4(|. Accounting for the in- 
version we find that, contrary to the conclusions of that 
paper, the results are most consistent with either deep 
minima or nodes close to the T-M direction. We predict 
the evolution of the C and n as a function of the field 
direction for different temperatures and fields. 

We consider a superconducting gap with the basis func- 
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FIG. 1: (Color online) Fermi surfaces and energy gaps in 
the two band model of iron-pnictides. The true structure is 
obtained by folding the FS sheets (blue) and SC gaps (red) 
along the sides of the square (full and dashed lines at M). 
The hole FS around T has an isotropic gap. The A\ g gap on 
the electronic FS may be isotropic, have deep minima along 
T-M line [panel (a), r = 0.45 in Eq. fl}] or nodes along the 
[100] as in panel (b), r = 1. We take the heat current along 
[100], as shown. 



tion in the A\ g representation, y(k x , k y ) = a + b(cosk x + 
cos k y ), where k x and k y belong to the unfolded Bril- 
louin Zone, Fig. HJa), and the ratio a/b determines the 
anisotropy of the gap. In that scheme a circular hole 
sheet is centered around T point, and the electron sheets 
are centered around M point and its equivalents. This 
choice of y{k x ,k y ) gives nearly isotropic order param- 
eter on the hole FS, A/,(0) = Ai, and, generally, an 
anisotropic gap on the electron FS, A e (0) = A 2 34 (</>), 

y e (4>) = (1 - r) =F r sin 20 , T for M,M' , (1) 

where <f> is the angle from the [100] direction, see Fig. [1] 
The gap has either minima along T-M line (r < 0.5, 
Fig. Q}i) or nodes close to T-M line for r > 0.5, which 
approach [100] and [010] directions for r = 1, Fig. \\]p. 
We consider these cases below. 

We compute the field-dependent specific heat and elec- 
tronic thermal conductivity by solving the equation for 
the quasiclassical Green's function g in the particle-hole 



(Nambu) space, as in Rcfs. [40|, |41[. To treat multi- 
ple bands we introduce <?n=i,2(R, k; e), with fc-parameter 
running over the hole FS (n = 1) twice (two equivalent 
hole bands), and over the two electronic FSs (n = 2) at 
points M and M'. In each band, given the Fermi veloc- 
ity v„(fc) and the vector potential A(R) of the magnetic 
field, we solve the transport equation for g n at energy e, 



e + -v n A(R) ) ? 3 - A„ - a n ,g n 



iv„ ■ V R g n = , 

~ (2) 

subject to the normalization condition c/ 2 = —ir 2 l. Here 
T3 is the Pauli matrix, and the spatial dependence of the 
order parameter is that of the Abrikosov vortex lattice, 



A n (R,k) = A n (k)^2C ly -^= 



e llyV ~ I x 
$ 



A 2 L 



A 



where &o{ z ) is the ground state wave function of a 
harmonic oscillator, x and y are in the plane normal 
to the field, and the magnetic length A 2 = hc/2eH. 
Two bands are mixed through the self-consistency on 
A„(R,fc) = TY, Sm ,k>, n ,V n , n ,(k,k')f n >(TL,k';e m ), where 
/„, the Gor'kov pairing amplitude, is the off-diagonal 
component of c/ n ; and on the impurity self-energy, a n , 
which is determined in the self-consistent ^-matrix ap- 
proximation for two bands [Hj], so that a(R;e) = 
nimpt s (R-', s). We take a negative inter-band pair hop- 
ping, V 12 (k,k') = -\V\y(k)y(k') which leads to the op- 
posite signs of gaps Ai and A2. 

We employ the extended Brandt-Pesch-Tewordt 
(BPT) approximation where the diagonal components of 
g n , the propagators g n and g n , are replaced by their spa- 
tial average, while the full dependence of /„ is ke pt. This 
approach is described and justified in Refs. (40l. 1-491. IHol] . 
and the results obtained from the self-consistent solution 
of the quasiclassical equations in a single band agree with 
those obtained using our method nearly perfectly. [HI 

Once the Green's function and the self-energies are de- 
termined, wc find the low-temperature specific heat from 



C(T,H) 
T 



^JL^osh- 2 — V JV n (e,T,H) 
T AT 2 2T ' 



(3) 

where A„(e,T,H) = -^(3mg n (e 7 k)) FS , is the DOS in 
each band and the angular brackets denote averaging over 
the corresponding FS. Similarly, k xx = K\ iXX + k 2 ,xx, 
where each Fermi surface contributes 



de e 2 , _ 2 £ 

x cosh — 

T 2T 2 2T 



(4) 



(vl x N n (T, H; k, e) T H . n (T, H; k, e)) 



FS 



and the transport scattering rate in each band is [41, 52 1 



2th,u 



(5) 

Here R indexes a retarded function, E„ = (d n )n, and 
V„ is the component of the Fermi velocity normal to H. 

Results. We take the e- and h- FSs to be cylinders of 
the same size, with |u n (fe)| = Vf. In the supplementary 
material we show that our conclusions are robust against 
modifications of the band structure, and only for the par- 
ticular case of high curvature of the electronic Fermi sur- 
face along the T-M line additional care is needed (53j . 
The unit for the magnetic field is Bo = fic/(2e£g), where 
Co = Vf/2irT c o is the in-plane coherence length, and T c o 
is the transition temperature in a pure sample. Highly 
anisotropic pairing states are affected by the disorder 
in each band. Hence below we present the results for 
the purely miro-band impurity scattering limit, with the 
normal state scattering rate, T/2ttT c q = 0.005, which 
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FIG. 2: (Color online) The anisotropy of the heat capacity 
(a,b) and thermal conductivity (c,d) for a gap with deep min- 
ima along F-M line (<j>o = 45°), as depicted in Fig. Ela). 
Except for the low-T, low-H regime, the minima in the gap 
are marked by the maximum of C. For heat transport, note 
the same trend, and almost complete absence of the two-fold 
anisotropy in k(4>o) at low H and T > 0.05T c q. 



FIG. 3: (Color online) The anisotropy of C/T (a,b) and k/T 
(c,d) for r = 0.55, gap with a pair of nodes close to the F-M 
line. The overall behavior is similar to that for deep min- 
ima, Fig. [2] The additional structure in C/T for near-nodal 
directions disappears already at low T, panels (a) and (b); 
Thermal transport shows additional structure to higher T, 
and the twofold component is dominant, panels (c) and (d). 



gives < 5% suppression of the transition temperature. 
We consider strong scatterers, phase shift 5 = n/2, and 
checked that smaller 6 and moderate inter-band scatter- 
ing do not perceptibly change our results. While the 
BPT method works well for nodal superconductors, its 
validity in systems with finite minimal gap A m m , is 
restricted to the regime where v//(2AA min ) > 1. |40j 
For r = 0.45 in Eq.jJT}, at the lowest field we consider, 
H = 0.02B , this ratio is about 2. Taking £ ~ 30l 
gives Bq ~ 35 T, so that the fields up to 14 T corre- 
spond to H < 0.4 Bo <C H C 2. The last inequality allows 
for non-self-consistent calculation of the order parameter 
suppression by the in-plane H, which is known to be an 



51|,|54| 



excellent approximation. [4C 

Figs. ©-(jU) show representative results for C and re 
as a function of the field direction at low fields for dif- 
ferent order parameter structures, r = 0.45,0.55, 1. The 
panels capture the qualitative behavior across the T-H 
phase diagram, with only quantitative changes as we go 
to higher fields and temperatures. The C and re profiles 
are slightly shifted vertically for clarity. 

Fig. [2] and Fig. [3] use the gap suggested by the major- 
ity of theoretical works, with cither minima [r = 0.45, 
Fig. HJa)] or closely spaced nodes [r = 0.55, inset in 
Fig- Eld)] in the T-M direction, <fio = 45° in our notation. 
The key feature is the inversion of the specific heat oscil- 
lations as the temperature is raised. At the lowest T, H 
the minima in C/T indeed occur for H along the T-M 
line, as expected from the semiclassical theory for zero- 
energy DOS, [13, [H, EH but this regime is narrow, and at 
higher T and H it is the maxima of the C-pattern that 
denote the minima/nodal directions. At low fields the 
inversion occurs at 0.05 < T/T c q < 0.1 for deep minima, 
Fig. HJa), and at even lower T/T ca < 0.05 for the nodal 
case, Fig. G2[a). The same inversion appears for all T at 



fields H > 0.W a , Fig. [^b) and Fig. E[b). In Ref. |4£g 
the C/T anisotropy was measured at T/T c ps 0.2. At this 
temperature the minima in C/T, observed at 4> = 0,90 
relative to the crystallographic axes, indicate deep min- 
ima at (f>o = 45° or nodes close to this direction as evident 
from comparison with the upper curves in Fig. [U[a,b) or 

Fig.EKb). 

In contrast, the experimentally observed pattern is 
not consistent with the nodes along [100] and [010]. 
Fig. IHa),(b) show minima in the C/T pattern for the 
field in the nodal direction, 4>$ = 0, only at T/T c q < 0.10; 
at higher T additional structure develops, followed by 
the inversion and the shift of the minima of C(0o) to 
<f>o =45°. This is not what was found in Ref. [46j ]. 

Hence possible gap structures are (a) minima along 
the T-M line, or (b) nodes close to this direction. Dis- 
tinguishing between the two by methods sensitive to the 
amplitude but not the phase of the gap is not straightfor- 
ward. While Fig. [3] clearly shows additional features (ab- 
sent in Fig. [2]) at the angles where |A e | has nodes, these 
features are washed out with increased temperature. For 
comparison, k/T shows nodal features at higher T (pan- 
els c,d) than C /T (panels a,b) . With increased scattering 
this structure smears out and largely vanishes when the 
nodes are lifted by disorder [HBI ] . 

Commonly, k(4>o) is decomposed into a twofold term 
due to difference in transport normal and parallel to the 
vortices [56| . and fourfold component due to gap struc- 
ture 42j, k = kq + K2 cos2</> + re4COs40. Both K2 and 
change sign in the T-H plane (ZJ, (Hif . The magnitude of 
«2 is similar in both Fig. [2] and Fig- [3J however, the four- 
fold K4 is much greater for the gap minima than for either 
nodal scenario. Observation of a large ratio k\/k,2 > 1 
would indicate minima and not nodes in the gap. 

Additional information is needed to distinguish be- 
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FIG. 4: (Color online) Anisotropy of C and k for r = 1, nodes 
along [100], as in 

Fig.Htb). The nodal directions 4> = 0,90°, 
are marked by minima of C only at T < 0.1T C . At higher 
T, H the angle-dependence shows a maximum of C for those 
directions, as well as additional structure. 



tween shallow and deep minima. In Fig. [2] for r = 0.45, 
(A m i n / A max = 0.1) the anisotropy in C is ~ 5% for our 
choice of the FSs. Setting r = 0.3 (A min /A max = 0.4) we 
found that this anisotropy drops below 1%, and its inver- 
sion occurs at 0.15 < T/T c q < 0.2 for similar fields. Also, 
above the inversion, the ratio K4/K2 is much smaller than 
that for r = 0.45. Of course, larger minimal gap should 
be evident from measurements even at zero field. 

Conclusions. We developed a framework for the cal- 
culation of the anisotropy in the heat capacity and ther- 
mal conductivity of two-band superconductors under ro- 
tated magnetic field, presented the results for several 
models of pnictides with anisotropic Ai g ('extended s') 
gap symmetry, and compared them with an experiment 
on the Fe(Se,Te) system. We identify either minima 
(A min /A max < 0.4) or nodes in the gap on the electronic 
FS along the T-M line, contrary to Ref. [Ifjj]. We predict 
that comparison of the fourfold and twofold term in the 
anisotropy of thermal transport will help distinguish be- 
tween the two scenarios. Experiments in a wider T-H 
range along with the calculations based on realistic FS 
arc clearly forthcoming. Our work lays the foundation for 
determining the gap structure from these measurements. 

We thank A. V. Chubukov, P. J. Hirschfeld, Y. Mat- 
suda and T. Shibauchi for valuable discussions, and ac- 
knowledge partial support from DOE Grant DE-FG02- 
08ER46492 (I. V.) and Aspen Center for Physics. 

Note added: After our manuscript was submitted sim- 
ilar conclusions were reached independently in Ref. (H3] 
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Supplement to "Nodes vs. minima in the 
energy gap of iron-pnictides from field-induced 
anisotropy" 

We show that our main results presented in the 
Letter are robust against the changes in the Fermi 
surface shape and the changes in relative magnitude 
of the gaps on different Fermi surface sheets. We 
present the anisotropy of the heat capacity, C(4>o), for 
the elliptical Fermi surfaces and for pairing interac- 
tions resulting in a significant difference in the gap 
magnitudes on the hole and the electron Fermi sur- 
face sheets. 

In our Letter we considered a simple model for the 
Fermi surface (FS) of the iron-based superconductors 
that consists of cylindrical Fermi surfaces both at the 
center of the Brilloin Zone (hole-like) and at the M and 
M' points (electron- like). While this is likely an adequate 
approximation to the electronic structure of the Fe(ScTe) 
compound where the only experimental data for the spe- 
cific heat anisotropy are available (4(§, broader theoret- 
ical questions include understanding for what classes of 
non-circular Fermi surfaces observed in the members of 
the iron superconductor family 0] the results remain fun- 
damentally the same, and under what circumstances we 
can expect significant deviations from the behavior out- 
lined in our Letter. The main point of interest is the 
sign of the field-induced anisotropic component of the 
specific heat in the T-H plane. An additional focus is 
on the detailed angle profile of the specific heat for dif- 
ferent Fermi surface shapes. The results presented below 
strongly suggest that, in most situations, our main con- 
clusions regarding the determination of the gap structure 
remain valid, although the magnitude of the oscillations 
and detailed angle-dependence arc sensitive to the varia- 
tions of the Fermi surface. We explicitly identify the cases 
where we can expect the results to be substantially mod- 
ified compared to our model predictions and be strongly 
dependent on the band structure for a specific material. 

There are three effects that we need to address. The 
first is the weak isotropic corrugation of the Fermi sur- 
face along the c-axis. In our previous work for a single 
band system with vertical line nodes we showed that such 
corrugation has a minimal effect on the location of the 
inversion line Q. The same holds true for the current 
model with either nodes or minima. More complex fea- 
tures may be expected if there is substantial modulation 
of the gap along the c direction: This is conjectured to 
happen in some of the compounds [H, 0] , but at present 
there is no universal agreement on what form the gap 
may take, and we postpone the analysis of this issue un- 
til a more detailed experimental and theoretical analysis 
narrows the possible gap structures. In the absence of 
a strong gap variation along the c direction our results 
hold, and we can, without loss of generality, consider a 



two-dimensional model. 

The second aspect that can potentially affect the vari- 
ation of the specific heat with angle is the difference in 
the gap magnitude on the electron and the hole sheets 
of the Fermi surface. Experiments indicate that one hole 
sheet may have a gap with the magnitude only half that 
on the other hole and electron Fermi surface sheets 0, 0] ■ 
We model an extreme scenario when both hole pockets 
have a reduced gap. Fig. [5] shows that our results are 
very weakly dependent on the ratio of the maximal gaps: 
only around T/T c < 0.05 and at moderately high fields 
there is a difference in the profiles. Over most of the pa- 
rameter range it is only the anisotropy of the gap that 
matters, and therefore changing the overall scale factor 
leads to a slight change in the number of available quasi- 
particle states, but not to a change in the pattern of the 
dependence on the direction of the field. 

The third, and potentially the most important effect is 
that of the electron Fermi surface shape. When the Fermi 
surface is not circular, it is the combined anisotropy of 
the Fermi surface and the gap that determines the angle- 
dependence of the specific heat and the thermal conduc- 
tivity. For the density of states at zero energy the Fermi 
surface effects have been investigated for borocarbide su- 
perconductors and for a model single Fermi sur- 
face at finite temperatures Fig. [6] shows the results 
for the case of pnictide-like two-band systems with the 
electron Fermi surfaces in the shape of an ellipse with 
the aspect ratio of 2, and the gap with deep minima, 
r = 0.45 in A e ((f>) = A 2 [(l - r) - r sin 20]. This is the 
simplest and most generic implementation of the Fermi 
surface anisotropy that does not rely on the material- 
specific calculations. Since the gap on the hole sheets is 
nearly isotropic the detailed shape of that Fermi surface 
does not affect our results. 



A, /T , 




FIG. 5: (Color online) Heat capacity anisotropy for differ- 
ent amplitudes of the gaps on hole and electron FSs. Pan- 
els (a,b): low-T, zexo-H gap amplitudes Ah = 1.57 T c q, 
A e = — 2.48T c o. Panel (c,d): much greater gap amplitude 
on the electron FS: A h = 1.25 T c0 , A e = -3.12 T c0 . 
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T/T =0.10 

_^ r d2, , , 



H=0.02 




H=0.50 



FIG. 6: (Color online) The anisotropy of the heat capacity 
for an elliptical electronic FS with the aspect ratio of 2 and 
the gap with deep minima for different fields (upper panels 
H — 0.02Bo, lower panels H = 0.5£>o) and temperatures 
from T/Tco = 0.02 to T/T c0 = 0.15. The contribution from 
each Fermi surface (upper and lower solid lines in each panel) 
depends on whether the minima in the gap are on the flat or 
strongly curved part of the ellipse. Dashed line is the average 
of the two Fermi surfaces. Dotted lines are the results of 
calculations for a circular FS. 



With finite ellipticity the two previously degenerate 
Fermi surfaces at point M become distinct as shown in 
the rightmost column, so that the minima of the gap 
located along the T-M line now occur along one of the two 
principal axes of the ellipse, see Fig. [5] The two scenarios 
correspond to the gap minima at the sharp end of the 
ellipse (the outer, "flower" -like, FS in the folded zone, 
shown by upper curves in Fig. |6]), and the gap minima 
along the flatter part of the Fermi surface (inner electron 
FS sheet, lower curves in Fig. [6]). In the latter case, the 
resulting profile of the variation of the heat capacity with 
the angle strongly resembles that for the cylindrical FS. 
Some additional structure is seen at temperatures close 
to the inversion line, and the inversion itself at low T is 
pushed to slightly higher fields. 

In contrast, when the minima/nodes are in the regions 
of the Fermi surface with high curvature, the inversion 
of the angle-dependent oscillation pattern occurs at very 
low values of T and H, in agreement with Ref. Left 
panel of Fig. [6] demonstrates that this inversion already 
occurred even at low field for T = 0.02T C . In this case 
small changes of Fermi surface curvature around the loci 
of the gap minima do affect the angle-dependence, and 
therefore for systems where this curvature is large at the 
intersection with T-M line, a detailed material-specific 
calculation may be needed. 

Thermal conductivity anisotropy exhibits similar 
changes. For the gap with deep minima at the flat parts 



of the ellipse, the amplitude of the fourfold term rela- 
tive to the twofold component is large; while for the gap 
minima at FS regions with high curvature the twofold 
component dominates. 

In some of the pnictide materials, such as 122 As-P se- 
ries, the electronic Fermi surface is closely approximated 
by the ellipse in the basal plane that undergoes 90° ro- 
tation as a function of the c-axis wave vector (l2l |. In 
this situation the measured anisotropy is essentially the 
average of the contributions from the two possibilities 
discussed above. This is also plotted in Fig. [51 and direct 
comparison with the cylindrical case shows overall simi- 
larity. There is a more complex angle-dependent struc- 
ture at some temperatures and fields, and at very low 
temperatures (T/T c < 0.05) the higher inversion field 
gives a narrow range where the oscillatory contributions 
have opposite sign for the circular and averaged ellipti- 
cal cross-section, see lower left panel. However, by the 
experimentally relevant temperature of T f=a 0.2T C the 
inversion already occurred for both cases. Therefore our 
main finding that the measurements at and around this 
temperature indicate the minima or the nodes along the 
T-M direction remains robust in this case. 

Importantly, we checked that for the elliptical electron 
Fermi surfaces and nodal gap the profiles of the thermal 
conductivity and the specific heat as a function of the 
field angle exhibit complex features qualitatively similar 
to those discussed in our Letter and not seen in experi- 
ment. 

Our conclusion therefore is that the main features of 
our results presented in our Letter remain valid in most 
situations when changes in the Fermi surface shape, but 
not its topology, are taken into account. Weak three- 
dimensionality or the Fermi surface and the variation in 
the relative gap amplitude between the electron and the 
hole sheets essentially do not affect the angle-dependence 
of the specific heat. The change in the Fermi surface 
shape has a greater effect, and detailed analysis is needed 
for the materials where the electron Fermi surface has 
high curvature at the intersection with T-M line for most 
values of the c-axis wave vector. Without this caveat we 
find, however, that within a wide range of parameters the 
inversion of the net measured anisotropy of the specific 
heat still occurs below T = 0.2T C , and therefore the ex- 
perimental results of Ref. [l| are indicative of the minima 
or the nodes of the gap along the T-M line. 

Of course, in some members of the large family of 
iron-based superconductors the Fermi surface may have 
not just a different curvature but an altogether differ- 
ent topology: this may happen with severe overdoping 
when one of the sheets disappears altogether [13[ or 
when additional order leads to a reconstruction of the 
Fermi surface as suggested, for example, in Ref. [14J. In 



these specific cases a separate investigation of the result- 
ing pattern is probably warranted, but our results are 
generally applicable to many compounds where the Fermi 
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surface consists of the hole sheets near the T and (perhaps 
anisotropic) electron sheets near the M (or X) points in 
the Brillouin Zone. 
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